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I ^ ' Abstract. Pentaquark states in lattice QCD probably lie close in energy to two particle 

scattering states. Correctly identifying the resonant state is a challenging, yet tractable, problem 
given the terascale computing facilities available today. We summarize the initial round of 
exploratory lattice calculations and discuss what should be accomplished in the next round. 



o 



X 



1. Introduction 

The experimental evidence both for and against the existence of a narrow S=l baryonic 
resonance around 1540 MeV is summarized elsewhere in these proceedings 0l2l- If confirmed, 
' this resonance lies 5-10% above the KN scattering threshold with an allowed strong decay. The 

implication for Lattice QCD theorists is that proper identification of this resonance in Lattice 
QCD simulations (even in the quenched approximation) should be substantially more difficult 
than the standard calculation of the low-lying mass spectrum of hadrons stable against strong 
decays. Fortunately, techniques are being developed to study hadronic excited states, scattering 
■ states and decays on the lattice. That this narrow resonance lies so close to the scattering 

threshold may make this one of the easier scattering problems to study. 

Initial Lattice QCD simulations studied the 5=1 ground state energies for various choices of 
total angular momentum J, isospin I and parity P IS El El 13 El El El • From tabled we see 
I that the lowest angular momentum (J=^) isosinglet and isovector channels have been studied, 

but not isotensor (1=2). In all these studies, the ground state energies were generally consistent 
with a KN threshold scattering state. If the first excited state energies were extracted and did 
not seem consistent with the next higher scattering state, this indicated evidence for a possible 
pentaquark state. 

The pentaquark operators used to create these quantum numbers on the lattice generally had 
all five quarks at a single spatial point. In this case, the spatial symmetry enables a continuum- 
like operator construction with a reasonable expectation of correctly identifying the of the 
operator. Table ^ indicates various strategies for constructing the operators. 

2. Hadronic scattering states on the lattice 

Hadronic correlation functions in Lattice QCD behave as a sum of decaying exponentials 

oo ^ 

c,(t) = {om o\m) = E ^ exp (-E,,t) ^ (1) 

where the operators Oi and Oj have the same quantum numbers (/, Iz, J, Jz, P, S, ■ ■ ■). Since the 
statistical noise grows exponentially in t it is very difficult to reliably estimate the contributions 
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Table 1. Current summary of lattice pentaquark 
spectroscopy (adapted from (llj. 
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Figure 1. Volume dependence of 
KN scattering. 



of even the first excited state {n=\) unless its energy Ei is not much larger than the ground 
state Eq and a judicious choice of operators Oi can be found which have much better overlap 
with the excited state {\Zii\ S> |-^jo|). This has been the basic approach of most of the studies 
in table ^ 

Assuming that an excited state energy can be reliably estimated, a resonance can be 
distinguished from a scattering state by varying the spatial volume. In a periodic box of size L^, 
momentum components are quantized: pi = riil'n jL. By varying L, energies of isolated moving 
hadronic resonances vary predictably: E']^{p) = E^i^i) ^ {nlix j Lf . The energy of a two-particle 
scattering state will also vary according to the scattering length cq: 
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EK{p)EN{q)L^ 



2.834^+ 6.375 



(2) 



Figure n illustrates the expected volume dependence of KN scattering for various isospin I and 
orbital angular momentum L channels where the scattering lengths have been experimentally 
determined. However, the energies of isolated resonances at rest do not vary in sufficiently large 
volumes where finite volume effects are not important. Thus, scattering states and resonances 
can be distinguished in principle, although a sufficiently wide range of volumes are seldom 
available in practice. 



3. Group-theoretic methods for hadron spectroscopy 

A better method for isolating hadronic excited states is a variational method where a large basis 
of A'^ different operators Oi are identified and an NxN Hermitian matrix of correlators C(t) 
is constructed. The matrix elements of C{t) are given by equation If we make a Schur 
decomposition of C(t) = ZA*Z^ and truncate the sum in equation to the first N terms, 
then we can make an ansatz that the N eigenvalues of A determine the N lowest energy levels: 
Xn = exp{—En)- In practice, it is better to find the solutions to the following equation, called 
a matrix pencil: 
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C(to) 



(3) 



The problem of isolating excited states is now reduced to finding a sufficiently large basis of 
operators such that some linear combination will overlap strongly with each state of interest. 

Over the past few years, we have developed a group-theoretic technique for constructing 
baryon operators that transform irreducibly under the lattice symmetry group jl2l 1131 ITU I15j . 
We have extended our technique to pentaquarks by constructing operators that are linear 
superpositions of gauge invariant terms of the form 

= c.-. (m)l {Df% (sr*);;, (ot^i^i (m)l w 

where A,B,C,D,E indices label quark flavor; a,b,c,d,e indices label color; a,/5, 7,/i,j^ are 

Dirac indices; tp {tp) indicates a smeared quark (antiquark) field; and Dj denotes the p-link 
covariant displacement operator in the j-th direction. The tensor Cabcde indicates one of three 
linearly independent color contractions needed to produce a color singlet. For the displacements 
considered here, the simplest choice Cabcde = ^abc^de was sufficient to produce a complete set of 
operators. 

For pentaquarks at rest, our operators transform as irreps of the double-covered octahedral 
group Oh- There are four two-dimensional irreps Gig,Giu,G2g,G2u and two four-dimensional 
irreps Hg and Hu. Continuum spin assignments J for lattice states must be deduced from 
degeneracy patterns across different Oh irreps jl6j . 

So far, we have constructed complete sets of two types of operators. The first type, called 
single-site, has all four quarks and the antiquark at a single site. The second type, called singly- 
displaced, has all four quarks at one site and the antiquark displaced by p sites along one of the 
six spatial directions. Table [21 shows the maximal number of linearly independent /z=0, 5=1 
pentaquark operators {uudds) for each isospin I and operator type. 

Table 3. Number of irreducible representa- 
tions for various /z=0, S=\ pentaquark op- 
Table 2. Number of linearly independent erators. Opposite parity irreps, e.g. Gig and 
lj;=0, S=l pentaquark operators. Giu occur with equal frequency. 

single-site singly-displaced 
Irrep 1=0 1=1 1=2 1=0 1=1 1=2 

Gi 19 32 12 76 138 54 
G2 3 9 3 44 92 36 
H 19 37 15 120 230 90 



Table El shows the maximal number of distinct irreducible representations of Oh that can be 
constructed for each operator type and isospin /. Irreps of opposite parity occur with equal 
frequency so the number of Gi irreps is the same for Gig and Giu, etc. Thus, if the observed 
S=l resonance at 1540 MeV is /=0 and J=^ then we can use up to 19+76=95 different Gi 
operators of either parity to construct our correlator matrices. However, if neither the single-site 
nor the singly-displaced type operators have a significant overlap with the resonance state, then 
even more operators will need to be constructed using different displacement patterns. 

The next step in our program is to extend our existing software which performs the needed 
Wick contractions to compute the elements of baryon correlation matrices to include pentaquark 
operators. We hope to report on initial quenched simulations in the near future. We expect 
these pentaquark resonances will have a relatively large spatial extent and heavy mass, so the 
use of improved actions on anisotropic lattices will be important. Once we have identified 



Operator Type 1=0 1=1 1=2 

single-site 240 460 180 

singly-displaced 1440 2760 1080 



linear combinations of operators which overlap strongly with the desired resonances, it may be 
worthwhile to use the same operators to compute correlation matrices on dynamical isotropic 
lattices as well. 

4. Conclusion 

The lattice studies completed to date have most likely seen the threshold KN scattering state 
in the 1=0,1 1^=^ S=\ channels. In some cases, there have been tantalizing hints of an excited 
state which have been interpreted as either scattering states or resonances. Until the lowest three 
or more excited states can be cleanly extracted over a range of volumes, it remains unlikely that 
lattice calculations can confirm or deny the existence of pentaquark resonances in QCD. The 
best method for extracting excited states is the calculation of large correlator matrices. The 
group-theoretic operator construction technique can be used to produce the necessary large basis 
of operators. 
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